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We present the results for the masses of the doubly heavy baryons Eqq/ and £Iq/q where Q, Q' = b, c obtained 
in the framework of the simple approximation within the nonperturbative string approach. 



Doubly heavy baryons are baryons that con- 
tain two heavy quarks, either cc, be or bb. Their exis- 
tence is a natural consequence of the quark model of 
hadrons, and it would be surprising if they did not ex- 
ist. In particular, data from the BaBar and Belle col- 
laborations at the SLAC and KEK B-factories would 
be good places to look for doubly charmed baryons. 
Recently the SELEX, the charm hadroproduction ex- 
periment at Fermilab, reported a narrow state at 
3519 ± 1 MeV decaying in A+K~ir + , consistent with 
the weak decay of the doubly charged baryon S+ Q . 
The candidate is 6.3er signal. Another candidate par- 
ticle is related by the replacement of a down by an up 
quark, so the mass difference is expected to be similar 
to that of the proton and neutron. SELEX, however, 
sees a mass difference 60 times larger ||). Whether 
or not the states that SELEX reports turn out to 
be the first observation of doubly charmed baryons, 
studying their properties is important for a full un- 
derstanding of the strong interaction between quarks. 

The purpose of this talk is to present the re- 
sults of the calculation || of the masses of the doubly- 
heavy baryons obtained in a simple approximation 
within the nonperturbative QCD (for a recent review 
see j| and references therein). The starting point of 
the approach is the Feynman-Schwinger representa- 
tion for the Green function of the three quarks propa- 
gating in the nonperturbatibe confining background. 
The role of the time parameter along the trajectory 
of each quark is played by the Fock-Schwinger proper 
time. The proper and real times for each quark re- 
lated via a new quantity that eventually plays the role 
of the dynamical quark mass. The final result is the 
derivation of the Effective Hamiltonian, see Eq. (|J) 
below. In contrast to the standard approach of the 
constituent quark model the dynamical masses mi 
are no longer free parameters. They are expressed 
in terms of the running masses defined at the 
appropriate scale of y, ~ 1 GeV from the condition of 
the minimum of the baryon mass Mg as function of 



dM B {mi) 
dm; 



= 0. 



(1) 



Technically, this has been done using the einbein 
(auxiliary fields) approach, which is proven to be 
rather accurate in various calculations for relativis- 
tic systems. 

This method was already applied to study 
baryon Regge trajectories S and very recently for 
computation of magnetic moments of light baryons 
0. The essential point of this talk is that it is very 
reasonable that the same method should also hold for 
hadrons containing heavy quarks. As in we take 
as the universal parameter the QCD string tension 
a. We also include the perturbative Coulomb inter- 
action with the frozen coupling a s (l GeV) = 0.4. 

From experimental point of view, a detailed 
discussion of the excited QQ'q states is probably 
premature. Therefore we consider the ground state 
baryons without radial and orbital excitations in 
which case tensor and spin-orbit forces do not con- 
tribute perturbatively. Then only the spin-spin in- 
teraction survives in the perturbative approximation. 
The EH has the following form 



H = 




(2) 



where Hq is the non-relativistic kinetic energy opera- 
tor, are the current quark masses and rrij are the 
dynamical quark masses to be found from (Q), and V 
is the sum of the perturbative one gluon exchange po- 
tential V c and the string potential V^ring- The string 
potential has been calculated in || as the static en- 
ergy of the three heavy quarks: String (r l, 1*2, r 3) = 
C-Rmin, where i? m i n is the sum of the three distances 
|t"i| from the string junction point. 

The baryon wave function depends on the 
three-body Jacobi coordinates 



(3) 



(4) 



(i, j, fc cyclic), where /i^ and fiij t k are the appropriate 
reduced masses 



m.irrij 
rii + m , 



(mi + mj)m k 



m, + nij 



mk 



(5) 



and /i is an arbitrary parameter with the dimension 
of mass which drops off in the final expressions. 

In terms of the Jacobi coordinates the kinetic 
energy operator Hq is written as 



i / d 2 



a 2 
d\ 2 



5__d_ K 2 {Q) 

2/x {dR 2 ^ RdR + R 2 



where R is the six-dimensional hyper-radius 



R 2 



(6) 



(7) 



and K 2 (fl) is angular momentum operator whose 
eigen functions (the hyperspherical harmonics) are 



K 2 (n)Y [K] =-K(K + 4)Y [K] , 



(8) 



with K being the grand orbital momentum. In terms 
of Y[k] the wave function ip(p, A) can be written in a 
symbolical shorthand as 

K 

In the hyper radial approximation which we 
shall use below K = and ip — ip(R). Note that the 
centrifugal potential in the Schrodingcr equation for 
the radial function x(i?) = R b / 2 ipK{R) with a given 
K 

(K + 2) 2 - 1/4 
R 2 

is not zero even for K = 0. 

In terms of the p and A the potential 

Vstring( 

ri,r2,r 3 ) has rather complicated structure. 
Let 9ijk be the angle between the line from quark 
i to quark j and that from quark j to quark k. If 
9ijk are all smaller than 120°, then the equilibrium 
junction position Y is 



ap + /3X, 



(9) 



where 



Table 1 

The constituent quark masses m.j and the ground 
state eigenvalues E (in units of GeV) for the var- 
ious baryon states. 



baryon 


mi 


m 2 


m 3 


E 


(qqq) 


0.372 


0.372 


0.372 


1.426 


(qqs) 


0.377 


0.377 


0.415 


1.398 


(qss) 


0.381 


0.420 


0.420 


1.370 


(sss) 


0.424 


0.424 


0.424 


1.343 


(qqc) 


0.424 


0.424 


1.464 


1.171 


(qsc) 


0.427 


0.465 


1.467 


1.146 


(ssc) 


0.468 


0.468 


1.469 


1.121 


(qqb) 


0.446 


0.446 


4.819 


1.085 


(qsb) 


0.448 


0.487 


4.820 


1.059 


(ssb) 


0.490 


0.490 


4.821 


1.033 


(qcc) 


0.459 


1.498 


1.498 


0.904 


(sec) 


0.499 


1.499 


1.499 


0.881 


(qcb) 


0.477 


1.524 


4.834 


0.783 


(scb) 


0.517 


1.525 


4.834 


0.759 


(qbb) 


0.495 


4.854 


4.854 


0.593 


(sbb) 


0.534 


4.855 


4.855 


0.570 



rrn + rrij 



1 4£ + (3-t 2 )cotx 
73 1+7 2 
p Z-t 2 



2>jjLij 2A sin x I + t 2 ' 



where 



2Acosx ■ 



3// 



ij,k 



and x is the angle between p and A. If 9ijk is equal to 
or greater than 120°, the lowest energy configuration 
has the junction at the position of quark j. It can be 
easily seen that dependence on m, in equations (||) is 
apparent and Y does not depend on quark masses as 
it should be. 

In what follows the string junction point is 
chosen as coinciding with the center-of-mass coor- 
dinate. Accuracy of this approximation that greatly 
simplifies the calculations was discussed in [9. Av- 
eraging the interaction V — V c + \4tring over the six- 
dimensional sphere one obtains the Schrodinger equa- 



tion for x(R) 

d 2 x(R) 
dR 2 



2(i 



R 8^R 2 



X (R) = 0, 
(10) 



where Eq is the ground state eigenvalue and 
2a s 16 1 / , 



(ii) 



Table 2 

Masses of doubly heavy baryons 



State 


this work 


lit] 


[HI 


[12 


[13] 


E{qcc} 

n{scc\ 


3.64 
3.82 


3.70 
3.80 


3.71 
3.76 


3.66 
3.74 


3.48 
3.58 


E{qcb} 
n{scb} 


6.93 
7.10 


6.99 
7.07 


6.95 
7.05 


7.04 
7.09 


6.82 
6.92 


E{qbb} 
Q{sbb} 


10.14 
10.31 


10.24 
10.30 


10.23 
10.32 


10.24 
10.37 


10.09 
10.19 



15tt 



We use a 



f2 



(12) 



0.009 GeV, m 

M 



(0) 



= 0.15 GeV 
0.17 GeV, to 



Q.s 
(0) 



0.39, m£ ] = 
1.4 GeV, and 



m^> = 4.8 GeV, slightly different from M. We solve 



Eq. (10) by the variational method introducing a 
simple variational Ansatz 



X (R) ~ R 5/2 e-» p2R2 



(13) 



where p is the variational parameter. Then the three- 
quark Hamiltonian admits explicit solutions for the 
energy and the ground state eigenfunction: E ss 

min E(p), where 
p 

E(p) = (x\H\x)=3p 2 

3 R 
IV 2 



15 



(14) 



_6_ 
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We first solve Eq. (0) for the dynamical 
masses mi retaining only the string potential in the 
effective Hamiltonian (|^). This procedure is in agree- 
ment with the strategy adopted in Ref . H . Then we 
add the perturbative Coulomb potential and solve 
Eq. (||) to obtain the ground state eigenvalues Eq. 
The results for various baryons are given in Ta- 
ble 1. The dynamical values of light quark mass 
m q ~ \fo — 400 - 500 MeV (q = u, d, s) qualita- 
tively agree with the results of Ref. obtained from 
the analysis of the heavy-light ground state mesons. 
For the heavy quarks (Q = c and b) the variation in 
the values of their dynamical masses toq is marginal. 
This is illustrated by the simple analytical results for 
Qud baryons. These results were obtained from the 
approximate solution of equation 



dE 
dp 







(15) 



in the form of expansion in the small parameters 

f = and ol s . (16) 

m Q 



Omitting the intermediate steps one has 

1/4 



Eu = -n n [ - 



\ + A-£--B-a a 



(17) 



1/4 



m i = V^\-) (l-A-Z + B-a s + ...), (18) 



m Q = mg } (1 + 0(e, a 2 a , a s £) + ...), (19) 
where for the Gaussian variational Ansatz ( |l3| ) 

A = ^Z1(») 1 '\ 0.24. B-i±*J*~ 0.42. 
2 \n J 18 V 7T 

(20) 

Note that the corrections of the first order in £ and a s 
are absent in the expression ( |l9| ) for toq. Accuracy 
of this approximation is illustrated in Table 1 of Ref. 

§ - 

To calculate hadron masses we, as in Ref. || , 
first renormalize the string potential: 



string 



string 



C, 



(21) 



where the constants Cj take into account the residual 
self-energy (RSE) of quarks j|. In what follows we 
adjust the RSE constants Ci to reproduce the center- 
of-gravity for baryons with a given flavor. As a result 
we obtain C q = 0.30 GeV, C s = 0.15 GeV, C c ~ 
C b ~ 0. 

We keep these parameters fixed to calcu- 
late the masses given in Table 2, namely the spin- 
averaged masses (computed without the spin-spin 
term) of the lowest double heavy baryons. In this Ta- 
ble we also compare our predictions with the results 
obtained using the additive non-relativistic quark 



Figure 1. Mass of as a function of the running 
c-quark mass for a = 0.15 GeV 2 and a — 0.17 GeV 2 . 
The masses are given in GeV. Bold points refer to the 
case mi ' = 1.4 GeV 
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model with the power-law potential |10|, relativis- 
tic quasipotential quark model pT| , the Feynman- 
Hellmann theorem |l2| ] and with the predictions ob- 
tained in the approximation of double heavy diquark 
[ p"3[ . The change of a to 0.17 GeV 2 increases the 
mass of S+ by ~ 30 MeV. The perturbative spin- 
spin interaction introduces an additional shift of the 
E!+ mass ~ —20 MeV. Note that the mass of H+ is 
rather sensitive to the value of the running c-quark 
mass ra°, see Fig.l. 

In conclusion, we have employed the general 
formalism for the baryons, which is based on nonper- 
turbative QCD and where the only inputs are a, a s 
and two additive constants, C q and C s , the residual 
self-energies of the light quarks. Using this formalism 
we have also performed the calculations of the spin- 
averaged masses of baryons with two heavy quarks. 
One can see from Table 2 that our predictions are 
especially close to those obtained in Ref. (hJ using 
a variant of the power-law potential adjusted to fit 
ground state baryons. 
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